ON THE DECOUPLING OF LAYERED SUPERCONDUCTING 
FILMS IN PARALLEL MAGNETIC FIELD 



J. P. Rodriguez 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 and 
Dept. of Physics and Astronomy, California State University, Los Angeles, CA 90032. 



Abstract 

The issue of the decouphng of extreme type-II superconducting thin films (A^ ^ oo) with 
weakly Josephson-coupled layers in magnetic field parallel to the layers is considered via 
the corresponding frustrated XY model used to describe the mixed phase in the critical 
regime. For the general case of arbitrary field orientations such that the perpendicular 
magnetic field component is larger than the decoupling cross-over scale characteristic of 
layered superconductors, we obtain independent parallel and perpendicular vortex lattices. 
Specializing to the double-layer case, we compute the parallel lower-critical field with 
entropic effects included, and find that it vanishes exponentially as temperature approaches 
the layer decoupling transition in zero-field. The parallel reversible magnetization is also 
calculated in this case, where we find that it shows a cross-over phenomenon as a function 
of parallel field in the intermediate regime of the mixed phase in lieu of a true layer- 
decoupling transition. It is argued that such is the case for any finite number of layers, 
since the isolated double layer represents the weakest link. 
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I. Introduction 



The study of layered superconductivity has been reinvigorated by the discovery of the 
high-temperature oxide superconductors.^ In the case of the Bismuth (Bi) and ThaUium 
(Tl) based compounds, for example, the superconducting coherence length perpendicular 
to the conducting planes is less than the separation between layers, and London theory 
fails. ^ Hence a Lawrence-Doniach (LD) type of description in terms of weakly Josephson- 
coupled layers becomes necessary.^ The following question then naturally arises: does a 
layer decoupling transition, aided perhaps by the introduction of a parallel magnetic field,^ 
occur in such systems in addition to or in place of the conventional type-II superconducting 
ones marked by the lower and upper critical fields,^ Hci{T) and Hc2iT)? 

In the absence of external magnetic field parallel to the layers, a layer decoupling 
transition is indeed predicted to exist theoretically at a critical temperature that lies 
above that of the intra-layer Kosterlitz-Thouless (KT) transition temperature, T^.^"^ This 
result is based on studies of the XY model with weakly coupled layers, which accurately 
describes a layered superconductor in the absence of fluctuations of the magnetic field;^'^'^ 
e.g., in the intermediate regime of the mixed phase found in extreme type-II superconduc- 
tors, where it is appropriate to take the limit — > oo for the London penetration length 
Xl- Also, recent experiments on Bi-based high-temperature superconductors find evidence 
for a superconducting transition at T^, where the c-axis resistivity vanishes, that lies a 
few tenths of a degree Kelvin above a second superconducting transition temperature T"'', 
where the resistivity in the a6-planes vanishes. ^^'^^ Hence both theory and experiment 
find an extraordinary regime in temperature, Tc = T^^ < T < = T^, where Josephson 
tunneling between layers exists while the layers themselves are resistive and show no intra- 
layer phase coherence. (It has recently been argued in ref. 13 that in fact marks a sharp 
crossover for the bulk case.) Last, it is worth pointing out that Monte-Carlo simulations 
of the three-dimensional (3D) XY model obtain similar behavior in the presence of a large 
magnetic field perpendicular to the layers. 

The nature of layer decoupling in the presence of magnetic field parallel to the layers, 
on the other hand, is less well understood. Early work by Efetov suggests that supercon- 
ducting layers decouple at high parallel fields B\\ > -B|(0) ~ ^o/(fj, where $o denotes 
the fiux quantum, 7 = {rric/mabY^'^ is the mass anisotropy parameter, and d denotes the 
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separation between layers. The latter calculation is based on a high-temperature series 
expansion analysis of the LD model. A recent study of the same model by Korshunov and 
Larkin that employs a Coulomb gas representation, however, finds no such layer decoupling 
in the high-field limit for temperatures below the decoupling transition. In this scenario, 
therefore, the layers remain effectively Josephson coupled up to the parallel upper-critical 
field, H^^{T). 

In this paper, we shall also examine the decoupling of layered superconductors in 
parallel magnetic field, but in film geometries of thickness much less than the in-plane 
London penetration length, and at temperatures near the zero-field decoupling transition 
at T*.^~^ Specifically, we consider thin films of extreme type-II layered superconductors 
(Xl cxo) in the intermediate regime (ifci <^ -ff ^ -^02) of the mixed phase, which can 
be described by a frustrated XY model with a finite number of weakly coupled layers. 
By working with the Villain form of the latter, ^'^^ we obtain first that the thermodynamics 
factorizes into independent perpendicular and parallel parts in the presence of magnetic 
field at arbitrary orientation so long as the perpendicular component, H±, of the latter 
is larger than the Glazman-Koshelev decoupling cross-over scale,^^ ~ ^o/d^l^- The 
perpendicular thermodynamics is characterized by the melting of two-dimensional (2D) 
vortices-^^ that are decoupled from the parallel Josephson vortices, as well as from the 
perpendicular 2D vortices in adjacent layers. This is a result of the fact that well-formed 
vortex loops traversing a few or more layers within the film are absent in the present limit 
of weak inter-layer coupling.^ In particular, the parallel Josephson vortices are unable to 
make connections between perpendicular 2D vortices in the same or in adjacent layer if 
the nearest-neighbor spacings between these perpendicular 2D vortices is much less than 
the zero-temperature Josephson penetration length, Aj(0) ~ ^d. This situation occurs 
precisely for perpendicular fields that satisfy H± » B^, as stipulated above. Also, since 
we first take the limit of extreme type-II superconductivity, we then have the inequality 
B^ ^ ~ $o/-^i- This means that the former requirement guarantees that the distance 
between perpendicular vortices is within the London penetration length, which in turn 
guarantees that magnetic screening transverse to the perpendicular field component is 
negligible. The parallel thermodynamics, on the other hand, is described by a LD model 
in parallel magnetic field H\\, with a heavily renormalized anisotropy parameter 7(T) that 

3 



diverges exponentially as T approaches the decoupling transition temperature from 
below. The renormalization down of the inter-layer coupling is due physically to the 
excitation of vortex rings^ (fluxons) that lie in between consecutive layers. 

Second, we compute the line tension of a single Josephson vortex in the simplest case 
of an isolated weakly coupled double-layer XY-model, where we find that the parallel 
lower-critical field if),'^(T) of the double-layer vanishes exponentially as T approaches 
from below. This result agrees up to a numerical factor with a recent calculation by the 
author of the same that employs an alternate "frozen" superconductor description of the 
Meissner phase in layered superconductors.^ Related results have also been obtained by 
Browne and Horovitz in the setting of long Josephson junctions^^ and by Horovitz via a 
fermion analogy for layered superconductors in parallel magnetic field. ^ The former coinci- 
dence is not surprising since the present double-layer type-II superconductor is equivalent 
to a (dynamical) long Josephson junction at zero temperature/^"^^ wherein no current is 
allowed to pass between the junction. In fact, the above analysis proceeds by first consid- 
ering the length of the vortex as imaginary time. Semi-classical quantum corrections to the 
energy of the fundamental sine-Gordon soliton,^^"^^ which corresponds to the Josephson 
vortex in the double layer, are then computed. Entropic wandering of the vortex therefore 
translates into quantum fluctuations of the soliton. Note that the wandering of Josephson 
vortices is equivalent to the excitation of double-layer fluxons, which are again the physical 
origin of this phenomenon. 

We next consider a one-dimensional lattice of Josephson vortices in the double layer. 
After adapting certain elements from the analysis of long Josephson junctions in external 
fleld^^"^^ to the "semi-classical" analysis discussed above in the case of a periodic array 
of sine-Gordon solitons, we are able to compute the reversible magnetization as a function 
of parallel magnetic fleld. Notably, we obtain a cross-over fleld b\{T) ~ ^o/(fi^{T), 
beyond which the magnetization displays a B^^ tail characteristic of both long Josephson 
junctions^^ and of layered superconductors in high parallel field, generally (see Fig. 1 
and ref. 27). Unlike long Josephson junctions, however, this cross-over field is much larger 
than the lower-critical field. Thus we find no evidence for field dependence in the layer 
decoupling transition temperature, T*, within the present "semi-classical" approximation, 
which is in agreement with the results of Korshunov and Larkin.-*^* Last, we argue that 
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such is the case for any finite number of layers, since the isolated double layer represents 
the weakest link. 

The remainder of the paper is organized as follows: in part II we introduce the frus- 
trated XY model in the Villain form, from which we derive the renormalized LD model 
for T < T^. The double-layer case is the focal point of section III, where we compute the 
parallel lower-critical field and the parallel reversible magnetization from the above LD 
model, in addition to the compressibility modulus of the corresponding vortex array and 
the effective inter- vortex interaction potential in the dilute limit. We then apply these 
results to the phenomenology of layered type-II superconductors in section IV, as summa- 
rized by Fig. 2. Last, we assess the validity of the present "semi-classical" approximation, 
as well as discuss the general case of N layers, in section V. 

II. Frustrated XY Model 

The object now is to understand extreme type-II superconducting films composed of a 
finite number N of weakly coupled layers in the presence of external magnetic field. In the 
intermediate regime of the mixed phase, where the magnetic field satisfies Hd <^ H <^ Hc2, 
the London penetration length is in general much larger than the inter-vortex spacing. 
Following Li and Teitel,^^ magnetic screening effects arc then negligible, and we may 
describe the system by a uniformly frustrated layered XY model with an energy functional 
given by 

N 

ExY = ^11 E E E {1 - cos[A^0(f, I) - A^(f, 0]} 

1=1 f l^=x,y 
N-1 

+ J±J2 E{1 - ^ + 1) - - Mr, I)]}. (1) 

1=1 f 

Here (j){f, I) is the phase of the superconducting order-parameter on the layered structure, 
where r ranges over the square lattice with lattice constant a, and / is the index for the 
layers separated by a distance d. We presume that the lattice constant a is larger than 
the size of a typical Cooper pair. The magnetic fliix threading the plaquette at site (r, I) 
perpendicular to the /j, = x,y,z direction reads = ($o/27r) ^ e^i^^A;^^^, where 
$0 = hc/2e is the flux quantum. Also, A^(;i!)(r) = 0(r + fi) — (f){r) is the lattice difference 
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operator. The nearest-neighbor Josephson couphngs are related to the respective masses 
in Ginzburg-Landau theory by 



J|| ={h'^ /2m\\a^){nsa^d), 



(2a) 
(26) 



where Ug labels the superfluid density. Notice that Jy is independent of the lattice constant, 
a, as required by scale invariance in two dimensions. Hence the anisotropy parameter, 
7' = (J||/J_l)"^/^, of the XY model is related to that of the mass, 7 = (m_L/m||)-^/^, by 



7 =7-- 
a 



(3) 



Throughout, wc will presume weak coupling between layers, 7' ^ 1. 

To compute the corresponding partition function Z = j 'D(f){r)e~^^^ ^^^'^ , we now 
make the usual (low-temperature) Villain substitution for the exponential factors above; 
i.e., e-^(i-^°«^) ^ (27r/?)-i/2 E~=-oo e^^^e-^'/^/^. After integrating over the phase-field,^ 
we then obtain the following dual representation equivalent to A'"-layered compact quantum 
electrodynamics (PQED) in the strong-coupling limit (modulo a constant) 

1 r 1 N-i 



exp 



1=1 



(4) 



where n^{r) is an integer link-field on the layered lattice structure of points r = (^jOj 
with jjL = x,y,z and n = (fixjUy). Here, P\\^± = J\\^±/kBT. Notice that the conserved 
integer field is conjugate to the superfiuid current — Af^ of the XY model (1) in 
the continuum limit. To proceed further, let us now decompose the parallel field n into 
transverse and longitudinal parts n{f,l) = n'{f,l) — n-{f,l) + n-{f,l — 1), where the 
transverse and longitudinal fields, n' and n_ , respectively satisfy the constraints V • n' = 
and V -fi- = Hz, with V = (A^;, Ay). We then take the customary potential representation 
n_ = — V$ for the longitudinal (inter-layer) field, which yields $(r, /) = YlfG^'^^f — 
f')nz{r', I), where 



G'(2)(r-) = 



B2 (27r)^ 4 — 2cos(A;2;a) — 2cos(A;ya) 
6 



(5) 



is (formally) the Greens function for the square lattice. In the limit of weak coupling, 
7' — > oo, the interlayer field n_ vanishes, which implies that n' is indeed an integer field. 
After making a suitable (lattice) integration by parts of the energy functional in Eq. (4), 
we then obtain the factorization Z — Zccfl^^^Z^Q for the partition function in the limit 
of weakly coupled layers, where 

4g= I] n^.(5[V-n'|^.,i]exp 

{n'(f,0} 



(6) 



is the partition function for the 2D discrete gaussian model, with the in-plane vector 
potential A = (A^^Ay) presumed to be in the gauge V • ^ = 0, while the inter-layer 
Coulomb gas factor reads 

f 1 ^ 

ZcG = Yl ^^P1 -^EEE^-^^"'^ - i)-n,{r,l)]G^^\r-r')[n,{r',l-l)-n,{f',l)] 

{n^(f,0} I " '=1 r,r' 

N-1 ^ N-1 ^ 

- i E E^^(^"'' o^.i^", - ^ E E^'(^"'' o k (7) 

1=1 r 1=1 f J 

with the fields at the boundary layers set to nz{r,0) = = nz{f,N). This factorization 
is consistent with the original layered XY model (1) that consists of N independent 2D 
XY models in this limit. For the more relevant case of 7' large compared to one, but 
finite, we now make note of the fact that these XY layers remain effectively decoupled in 
the presence of perpendicular magnetic fields that are larger than the Glazman-Koshelev 
decoupling scale, -"^^ ~ $o/Aj(0), where 

Aj(0) = 7'a = 7d (8) 

is the Josephson penetration length. We therefore argue on a physical basis that the above 
factorization prevails in the presence of Josephson coupling as long as the perpendicular 
field satisfies Hj_ ^ in thin films of layered superconductors. 

In the absence of external magnetic field, A^{r) = 0, each layer (6) thus undergoes 
a KT transition at /c^Tc ~ f <^||) while the inter-layer links n^(f) undergo an inverted 2D 
Coulomb gas binding transition at A;bT* = 47rJ|| in the limit of weak inter-layer coupling, 
7' ^ 00.^ (It is understood that the limit of vanishing Josephson coupling is taken before 
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that of vanishing perpendicular field.) The latter high-temperature transition, which oc- 
curs well inside the normal phase of each individual layer, corresponds to the decoupling 
of layers mediated by the binding of oppositely (n^) charged vortex rings lying in between 
consecutive layers.^ Notice this implies that Josephson coupling between resistive layers 
exists in the temperature regime Tc < T < T*!^'^^'^^ For 7' large but finite, the form (7) 
of the layered Coulomb gas ensemble indicates that each set of consecutive double layers 
is dielectrically screened by itself as well as by the neighboring N — 2 such double layers. 
Hence, they each can be considered in isolation from their neighbors as long as one makes 
the replacement /Jy — > e^~^/?||, where eo — 1 is the polarization of an isolated double layer. 
Since the latter is directly proportional to the fugacity, = exp(— 2777'^), of the Coulomb 
gas (7) at [see Eq.(ll)],^ and since the decoupling transition temperature is then given 
by ksT^ = 47r J|| [1 + (A^ — l)(eo — 1)] in the limit of weak inter-layer coupling, we obtain 
an implicit linear dependence for the corrections to the value of with the fugacity y*. 
Such a linear dependence agrees with the standard renormalization-group fiows that cor- 
respond to the 2D Coulomb gas.^° As expected, the above formula for also indicates 
that the decoupling transition temperature increases without bound with the number of 
layers. In particular, the former linear increase crosses over to an exponetial increase at 
-^0 ~ (^0 — l)""*^ layers, which is far beyond the 2D-3D cross-over point expected to occur 
at Nc/o ~ 7' layers for the layered XY model. ^-"^ This is then consistent with the fact that 
the bulk layered XY model exhibits only a 3D supcrfluid transition at the bulk Tc. Note 
that the present factorization into parallel and perpendicular parts is unable to obtain 
corrections for the value of Tc in the case of large but finite anisotropy parameters 7'.^^ 
Consider now Eqs. (6) and (7) in the presence of a homogeneous magnetic induction, 

with the parallel component directed along the y-axis, and with B± ^ to insure the 
decoupling between 2D perpendicular vortices and parallel Josephson vortices explicit in 
the previous factorization of (4).^^ This decoupling becomes evident if we choose the gauge 
Aj: — 0, Ay = b±x, and = —b\\x, where b\\ and b± are the parallel (7) and perpendicular 
(6) magnetic flux densities. In particular, each layer independently experiences the per- 
pendicular component B± = H± of the magnetic induction,^^ which sets the intra-layer 
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vortex density to be ny = \H±\/^q. The fact that ny ^ Aj^ » A^^ in the present hmit of 
extreme type-II superconductivity insures that magnetic screening effects transverse to the 
perpendicular field component can be neglected. Each layer will then independently follow 
the 2D vortex lattice melting scenario/^ with a melting temperature Tm < Tc. Since the 
issue of 2D vortex lattice melting has already been discussed extensively in the literature 
with respect to the phenomenon of high-temperature superconductivity/ we shall end the 
present discussion here and focus our attention below on the thermodynamics connected 
with the parallel component to the magnetic induction. 

We now derive the renormalized LD theory^ ~^ mentioned in the introduction. It is 
useful first to make the following Hubbard-Stratonovich transformation of the Coulomb 
gas ensemble (7):^"^ 

ZcG=fveif,l) Yl exp|-4^J^(w)^_^^J^^(f,;)[^^(f,;)_^^(f,;_l)] 

{nzir,l)} I 1=1 r 1=1 r 

N-1 N-1 ^ 

E T.^z{r,i)A,{f,i) - ^ E E"'(^"''0 • (10) 

1=1 f 1=1 r J 

Here ^(f, I) represents a real scalar field that lives on each layer. Suppose now that we 
operate in the low-temperature regime, T < T*, where the layers are coupled.^"^ Then 
inter-layer charges in the Coulomb gas ensemble are screened, which means that global 
charge conservation is no longer enforced. Following the standard prescription,^^ we in- 
dependently sum over charge configurations at each site, with the restriction to values 
= 0, ±1. In the limit that the fugacity 

yo = exp(-l/2/3^) (11) 

is small, we then find that the original Coulomb gas ensemble (7) is equivalent to a renor- 
malized Lawrence-Doniach model Zcg = J 'Dd{f,l)exp{—Ei,j:,/kBT) with energy func- 
tional 

N N-1 

Elb = ^11 E E 2^^^)' + E Ei 1 - ^°«[^(^"'' ^ + 1) - ^(^"'' - Mr, 0]}- 

1=1 f 1=1 f 

(12) 
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At the decoupling transition in particular, we have that the fugacity (11) is given by 
= exp(— 2777'^). Hence, the anisotropy parameter is renormalized up to 

y^=( —i—] ' = (87r)-V2e-^" (13) 

at the decoupling transition in the present LD model. Since the charges physically 
represent vortex rings (fluxons) that lie in between consecutive layers,^ we conclude that 
such excitations are responsible for the renormalization (13) of the mass anisotropy near 
the decoupling transition. In closing, we remind the reader that the above derivation of 
model (12) is valid only for fugacities (11) that satisfy yo -C 1; i.e., for temperatures 
T > JjkB. 



III. Double Layer 

We shall now consider the parallel thermodynamics associated with the renormalized 
LD model (12) in the special case of a double layer (AT = 2), which is analytically tractable. 
This case is very similar to a long Josephson junction^^"^^ restricted to pass no current 
between the junction. Although the vortex lattice that results from the infinite-layer LD 
model in parallel field has far more structure than the simple vortex array corresponding 
to an isolated double layer, we believe that it is sufficient to study the latter with respect 
to the issue of layer decoupling in general, since it represents the weakest link. 

For the special case of two weakly coupled layers in the presence of a homogeneous par- 
allel magnetic induction directed along the j/-axis, the partition function corresponding 
to the renormalized LD model energy functional (12) is expressible as 

ZcG = J Ve{f)Vcf>{r')exp J -hp' J dyLsG[<t>] - J d\'^{Vef I, (14) 



where e{r) = 2-^/^[e{f, 1) + 9{f, 2)]. Here 
-^sg[0] = I dx 



(15) 



dy i 2\dx 

represents the "Lagrangian" for the sine-Gordon model in one space (x) and one imaginary 
time {y — it) dimension, with a bare Josephson penetration length 

Ao = a(/3||/4yo)'/', (16) 
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while the effective dimensionless Planck constant is 



^F = 2//?||. (17) 

Notice that we have taken the continuum limit of the LD model (12), as well as made the 
change of variable 0(r) = ^(r, 2) — 9{f,l) — Ag{r). The integration over 6 on the right of Eq. 
(14) results in a trivial gaussian factor. Below, we shall exploit the quantum mechanical 
analogy suggested above for the nontrivial sine-Gordon factor in order to compute the 
parallel lower-critical field^ and the parallel reversible magnetization of the double-layer 
system. 

A. Single Josephson Vortex 

We now set ourselves to the task of computing the parallel lower-critical field, H^}i, of 
the double layer system, which is in general related to the free energy per unit length of 
a single Josephson vortex, ey, by-*^ H^-^ = 47r£||/$o- Let us therefore consider the effective 
sine-Gordon model (15) in the presence of a single Josephson vortex aligned along the 
y-axis; i.e., the homogeneous magnetic flux is set to fey = 0, while the phase-difference 
field winds once along the a;-axis; dxd(f)/dx = 27r. In the absence of thermal (or 
"quantum") fluctuations, the vortex line tension is given by the Ginzburg-Landau energy 

4 = ^Lso[ct>o] = (18) 
of the "static" fundamental sine-Gordon soliton 

0o(a;,y) = 4tan-ie^/'^°, (19) 
which is a solution of the field equation 



+ Ao^sin0 = O (20) 



obtained by minimizing the corresponding "action" L^Lgcf'/'o]?^'^'^^ and which represents 
the single Josephson vortex. For temperatures near the decoupling transition at T*, how- 
ever, vortex wandering is critical,® and entropic (or "quantum mechanical") corrections to 
the vortex line energy (18) must be accounted for. 
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We shall now include the effects of thermal wandering in the double-layer Josephson 
vortex (19) by first Wick rotating the y coordinate to a real time-like coordinate, y = 
it. Second, we observe that the free energy per unit length of the Josephson vortex is 
equal to the product of ^ Jy with the "quantum mechanically" renormalized "mass" of the 
fundamental sine- Gordon soliton. To obtain the latter, we shall employ the "semi-classical" 
approximation^"^'^^ generally valid in the limit Hf ^ 0; i.e., at low temperature. In 
particular, consider small deviations = 00 -h 0i from the "static" vortex configuration 
(19). Then integration by parts yields that (15) is approximately 

Lsg[(I>]=Lsg[(I>o] + IJ ^^^</'^(|^-^+Ao 'cos 00^01 (21) 

to second order in the deviation. Hence in the presence of the soliton, we obtain a spectrum 
of harmonic oscillators of the form (l)i[x,t) — ip{x)e'^^^, where 

^2 \ 
■^ + Ao'cos(^o lV' = ^V (22) 

It is well known^^'^^ that the spectrum corresponding to (22) is composed of a zero-mode 

ipb{x) = sech (x/Aq), L>Jb = (23) 

that lies within the gap of the continuum 

V;fe(a;) = e^'^^f/c + zAq Hanh(a;/Ao)], uJk = {k^ + Aq ')'^'- (24) 

In the "semi-classical" approximation,^^ the "transition amplitude" over a period of "time" 
To is given by the product 2'sg[1] = exp(— zTo£|j/A;_BT)nfc2;fc, where 



^fe = X] exp 



n=0 



-iiOkTo I n -I- ^ 



2 



(25) 



gives the corresponding amplitude of each harmonic oscillator. After Wick rotating back 
to imaginary time Ly = zTq, we observe that only the n = terms above survive the 
limit of a long vortex, Lj^ — > oo. Yet the ratio of the partition functions in the presence 
of a single Josephson vortex to that in its absence is in general related to the vortex line 
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tension, £||, by Zsg[1]/-^sg[0] = exp(— Lj/£||/A;bT). In the "semi-classical" limit, therefore, 
we obtain 

-ii = 4 + ¥^fE-*-E-») (26) 

\ k q / 

for the line tension, ^'^'^^ where the k sum and the q sum above correspond respectively 
to traces of the zero-point energy in the presence and in the absence of the fundamental 
sine-Gordon soliton. In particular, the presumption of periodic boundary conditions along 
X implies the quantization conditions kL^ + S{k) = 27rn and qL^ = 27rn, where 

5{k) = 2 tan-\l/kAo) (27) 

is the phase shift of the continuum states (24), and where n is any integer. Properly 
counting these states then yields that the difference in brackets, J^ni'^k ~<^q)i in Eq. (26) 
is equal to^^'^^ —{2TT)~^j^^dk{duJk/dk)5{k) in the limit La; — > oo. After introducing a 
momentum cut-off, r^^, and integrating by parts, we obtain 

£|l=£j[l-(47r/3y)-iln(Ao/ro)]. (28) 

Notice that the correction above to the low-temperature line tension is first order in ^i?, 
and that it can be interpreted as a renormalization to the "mass", A,^^, of the sine- Gordon 
model (15). Since a renormalization group exists, we may now express the line tension as 

£||=4J||/Aj, (29) 

and then iterate Eq. (28), which yields AJ^ = Aq ^exp[— (47r/3||eo)~^ln(Aj/ro)], or 

l-(47r/3||eo)-']-' 

. (30) 



A 



J 




Above, we have replaced the right-hand side of Eq. (28) by the previous exponential and 
included the dielectric correction eo to the on-site Coulomb-gas (7) potential, ln(Aj/ro). 
Employing the standard renormalization group result, eo = 1 + 0[{(3\\ — /3*)^/^], for 
the dielectric constant of the 2D Coulomb gas at (inverted) temperatures f3\\ just above 
the (inverted) transition temperature /?* = l/47r, we thus obtain that the renormalized 
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Josephson penetration length Aj diverges exponentially as it approaches the decoupling 
transition like 

Aj/a = Cexp[D/(/3y-/3,)V2]. (31) 

Here, C and D are non-universal numerical constants. 

In conclusion, the parallel lower-critical field = 47r£||/$o vanishes exponentially- 
fast near the decoupling transition of the double layer following Eqs. (29) and (31). 
Horovitz has obtained this result employing a fermion analogy for layered superconductors 
in parallel magnetic field. ^ A similar dependence has also been proposed by Browne and 
Horovitz for the lower-critical field of long Josephson junctions. Combining Eq. (2a) 
with the identity ($o/Al)^ = (27r)^(^^/2m||)ns for the bulk {N oo) in-plane London 
penetration length Xl yields 



from which we obtain the useful expression 







2 



(33) 

XATY ^^^^ 



Given that C ~ 7' in Eq. (31), which is consistent with expression (8) for the Josephson 
penetration length of the original XY model (1) at zero temperature, then we have by 
Eq. (31) that Aj(0)/Aj(T) ~ exp[— iI>/(/3|| — /3^<)^/^] near the decoupling transition. This 
result agrees up to a numerical constant with a previous calculation by the author of the 
same quantity using an alternative "frozen" superconductor model for highly anisotropic 
extreme type-H superconductors in the Meissner phase.^ Notice then that expression (33) 
for the parallel line tension is essentially independent of the lattice constant a, as expected 
by the 2D scale invariance of the LD model (1) for T > T*. 

Before we go on to consider an array of Josephson vortices in the next section, a few 
remarks are called for. First, notice that the point at which the soliton "mass", 2£||/J||, 
vanishes coincides with the layer decoupling transition. Similar effects are found in the 
case of the nonlinear a model in two space and one time dimensions, which describes the 
quantum 2D antiferromagnet. In particular, the quantum mechanically renormalized en- 
ergy of the corresponding topological soliton called a skyrmion vanishes precisely at the 
zero-temperature quantum critical phase transition into the quantum disordered phase 
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characterized by a spin-gap. ° Second, notice that we have essentially recovered the stan- 
dard renormalization group results for the KT transition^° via the present semi-classical 
quantization of the sine-Gordon soliton energy in one space and one time dimension. Fi- 
nally, also observe that the entropic correction to the line-tension in Eq. (28) indicates that 
the number of microstates per unit length a of a Josephson vortex in thermal equilibrium 
is Xj/tq. Given that ro ~ a, then Aj can be naturally interpreted as the effective width 
of the Josephson vortex. 



B. Array of Josephson Vortices 



Consider now the case of a nonzero homogeneous magnetic induction aligned parallel 
to the y-axis of the double layer; i.e., b\\ ^ 0. Then it is easily seen from Eq. (15) that the 
superfluid portion, 

G,-Gn = -ksT In j Vc^e'^^' I '^^^^^ , (34) 

of the Gibbs free energy"^^ is minimized with respect to fey at 6|| — I^oo 94'/9x. In other 
words, the average winding per unit length in any configuration of the phase difference 
between the double-layers is set by the magnetic induction. In the particular case of 
the low-temperature "classical" configuration, we then have that the parallel magnetic 
induction is related to the lattice constant ao of the corresponding array of Josephson 
vortices by 

h = ^- (35) 
ao 

Clearly, we expect qualitative differences between the thermodynamics of the the low- 
field regime, ao ^ Aj, and of the high-field regime, ao <^ Aj. Yet does a decoupling phase 
transition separate the two regions, as has been suggested in the literature?^ Below we give 
evidence for the existence of only a crossover^^ on the basis of a "semi-classical" analysis 
of the reversible magnetization (see ref. 27) and of the elastic compression modulus of the 
vortex array. 

We now set ourselves to the task of computing the parallel reversible magnetization^^ 

Mil ^--^( I (36) 

15 



of the double layer in the intermediate regime of the mixed phase, W^-^^ <^ By <^ H^^2j 
where Hn = B\\. Consider first the lowest order Ginzburg-Landau contribution 



G^g -Gn = Ly-^LsG[(t>o] 



(37) 



to the Gibbs free energy in powers of hp, where (t)Q{x) represents the whirling pendulum 
solution of field Eq. (20) with spatial period oq; i.e.,^^~^^ 



#0 

dx 



-dn 



X - Xq 



K 



(38) 



where the parameter k lies in the interval between zero and unity, and is set by the period 
ao following 

ao = 2KqkK{k'^). (39) 

Above, dn(tt|K^) represents the appropriate Jacobian Elliptic function, while K{k^) rep- 
resents the complete Elliptic integral of the first kind."^^ Conservation of energy, Eq = 
2Aq^(k~^ — 1), in the corresponding pendulum system yields 

2n 



= 2 



Jo 



2\ dx 



En 



2^11- 



We obtain, therefore, that the zero-order Gibbs free energy density is equal to 



LxLy 



2 E(k^) 



+ 1- K~ 



2 2 ' 



(40) 



k2 Kin"^) 

where -E(k^) represents the complete Elliptic integral of the second kind.^^ Standard 
manipulations-^^ then yield that the reversible magnetization (36) is given by 



-47rMii = H. 



cl 



1 



4 kK{k^ 



(41) 



where H'^^ = 47r£|j/^»o is the parallel lower-critical in the Ginzburg-Landau theory approx- 
imation. Notice that H^-^ naturally sets the maximum value of the diamagnetic magneti- 
zation (41) at zero magnetic induction (k, = 1). In particular, at low magnetic inductions 
ao 3> Ao, we have by (39) that = 1 — 16e~"'°^^°. We then obtain the limiting behavior 



-47rMii ^ H 



cl 



1 + 4e ^0 + 1 - — — 



Ao 



2 ao 



(42) 
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for the reversible magnetization (41). Hence, the low- field magnetization extrapolates to 
zero at 

II 

^°"7r2Aod' ^^^^ 
which defines a bare crossover field. At high- fields oq <^ Aq, on the other hand, (39) 

dictates that n = 7r~-'^ao/Ao. We then obtain that the limiting behavior for the reversible 

magnetization (41) is given by 

_o II 

in such case. This implies a By tail at high fields By » Bq that is characteristic of long 
Josephson junctions and of layered superconductors in general. ^^'^^ A plot of result (41) 
spanning both the high-field and low-field limits is shown in Fig. 1. Notice that the bare 
crossover field Bq is much larger than the Ginzburg-Landau lower-critical field H^^ in the 
present double-layer extreme type-II superconductor [see Eq. (51)]. This is qualitatively 
different from the case of a long Josephson junction, where Bq ~ -ffji- 

In analogy with the previous analysis of a single Josephson vortex, let us now consider 
the effect of "semi-classical" corrections to the reversible magnetization in parallel field 
(41). We again have a spectrum (21) of harmonic oscillators 0i(a;, f) = 'ip{x)e^^^ that 
satisfy the linearized field equation (22), but with a periodic configuration for the phase 
difference set by^°'^^ 

(45) 



1 , I X — Xq 

COS -00 = — sn 



2 \ kAq 

where sn{u\K^) represents the appropriate Jacobian Elliptic function. To be more specific 
the spatial factors of each oscillator satisfy Lame's equation,^-*^ 



+ 2Ao ^sn^ 



t(; = cc;>. (46) 



To make contact with the previous discussion of a single Josephson vortex, let us now 
focus our attention on the (bare) low-field regime By <^ Bq, where the parameter k is 
exponentially close to unity, since ao 3> Aq. This allows us to approximate the potential 



terms in Lame's equation [— ^-j + V{x)]%l) — (a;^ — Aq ^)'0 by 



vi.)^- j: 2a„-w( "--;;"°° ). (47) 
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where each term above corresponds to the potential associated with a fundamental sine- 
Gordon soliton centered at xq + nao. In general, the band structure corresponding to 
Lame's equation (46) is composed of a continuum and a zero-mode band separated by 
a gap.^^''^^ A curious feature particular to each potential term in Eq. (47), however, is 
its transparency;^^ i.e., the continuum oscillators (24) of the fundamental sine-Gordon 
soliton have no reflected wave component. Therefore in the present (bare) low- field limit, 
the upper continuum band is essentially the same as that of a fundamental soliton (24). 
Repeating the renormalization group arguments made in the previous section for the line 
energy of a single Josephson vortex then indicates that the entropic correction due to the 
latter continuum band can be accounted for by simply replacing Aq (16) with Aj (31) in 
the original Ginzburg-Landau free energy of the vortex lattice; i.e.. 



Gs-Gn = ^ I dy I dx 



2 / \ 2 



dy I 2\ dx } 2 



(48) 



with the lattice constant of the vortex array (38) set by Eq. (35). In general, however, the 
effects of the zero- mode band must also be included in the present semi-classical analysis. 
The corresponding states are given by the tight-binding anzats |fco) = Ylm e^^°°'°'^\n) in the 
present (bare) low-field limit, where {x\n) = •0b (x — xq — nao) is the (normalized) bound 
state (23) of the fundamental soliton located at the n*^ well. The hopping matrix element 
is therefore -to = (n| - V(x)|?i + 1) = (w^-Aq ^)(n|n + l). But {n\n + l) = 27re""o/'^o 
is much less than unity, which yields to — AQ)e~"'°/^° . This means that the zero-mode 
band has a spectrum u>o{ko) = 2ty^|sin(|A;oao)| that is exponentially narrow. By (25), 
the zero-mode band results in an entropic pressure contribution to the Gibbs free energy 
density given by 

^0 = ^^. }:r'iko) = -^t,^ • (49) 

feo 

Hence, the magnetization (36) acquires a diamagnetic correction —dPo/dB\\ of order 
g-ao/2Ao^ which is negligibly small in the present (bare) low-field limit. Eq. (42) indi- 
cates, however, that the low-field correction to the initial linear increase of the parallel 
magnetization varies as 6"°^°/^° in the Ginzburg-Landau regime. Unlike the case of a 
single Josephson vortex, then no obvious renormalization group appears to exist for the 
above entropic pressure contribution. 
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In conclusion, double-layer extreme type-II superconductors in parallel magnetic field 
are described by the effective Ginzburg-Landau free energy (48) , along with the boundary 
condition (35), in the bare low field limit B\\ <ti Bq of the intermediate regime, H^^ -C 
B\\ <^ -f^c2- This means that the reversible magnetization is determined by the original 
Ginzburg-Landau theory analysis [Eqs. (36) - (44)], where the bare Josephson penetration 
length Ao is replaced by the renormalized length Aj throughout. In particular, the true 
parallel cross-over field (see Fig. 1) of the double layer is given by 

b! = AilL (50) 

instead of by Eq. (43). However, our inability to find a renormalization group for the 
entropic pressure contribution (49) to the parallel magnetization suggests that the present 
renormalized Ginzburg-Landau theory result for — 47rM|| serves only as a strong lower 
bound in the critical regime. We therefore find evidence for at best a crossover as a 
function of magnetic field below the bare scale B^ , and no evidence for a decoupling phase 
transition at fixed temperature. Last, it is easily shown after employing relation (32) that 

^" - (51) 



-"cl 

This of course indicates that the crossover field is much larger than the lower-critical field, 
which validates a posteriori the assumption (36) that i/y = in the intermediate regime 
of the mixed phase. It also illustrates the qualitative difference between a double-layer 
superconductor and a long Josephson junction, -"^^ where bI ~ H^^-^. 

We shall close this section by computing the compression modulus of the parallel array 
of Josephson vortices, as well as the interaction energy between widely spaced vortices. The 
local change in the elastic free-energy density due to a local fluctuation Sny in the vortex 
density is given by 5/sg = ^d'^ fsG/dn'^){5nv)'^ , where 



JSG = 



«2 K(k2) 



(52) 



is the Gibbs free-energy density (48) modulo the — l^*! term [see Eq. (40)], which can 
be considered as part of the magnetic field energy. In general, the number of Josephson 
vortices per unit length along the a;-axis is ny = clq^ = 6||/27r. By differentiating the 
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first term in Eq. (41) once with respect to we thus obtain d'^fsc/driy = (8J||/d)(l — 
k,'^)[K{k,'^)]^/E{k,'^). Now the displacement field u{x) of the vortex array along the x- 
axis is related to the density fluctuation by Sny = —aQ^{du/dx). Therefore employing 
previous identities [(9a), (32) and (35)], we find that the elastic energy is given by 5/sg = 
^cii{du/dx)'^ , where the compression modulus reads 



Cll = TT -—T—{l-K,). (53) 



In the low-field limit By <^ Bll we then have that the array is exponentially soft,^^ with 
Cll e~°'°/^'^. Similar softening of the Abrikosov vortex lattice occurs in conventional 
superconductors near the lower critical field, but with the Josephson penetration length 
A J replaced by the London penetration length A^-^^ On the other hand, the high-field 
limit ^ bI yields that cn = (47r)~^(d/Ai)^i?y , which is a formula characteristic of the 
elastic moduli in extreme type-II superconductors (Al ^ oo) generally.^ 

Last, we may define the interaction energy between two well separated vortices by 
subtracting the line tension (29) from the free energy (52) per unit length of a single 
period in the array of Josephson vortices; i.e., the repulsive interaction energy per unit 
length is given by 

v{ao) = daofsG - £|| (54) 
in the low-field limit ao ^ Aj, which after some analysis yields v(ao) = J\\Xj^(l — k^), or 

vix) = 16J||A7'e-l^l/^^. (55) 

Notice that the Josephson penetration length Aj acts as the screening length for the 
interaction between Josephson vortices instead of the London penetration length, which 
plays the same role in conventional type-II superconductors. 



IV. Phenomenology 

We shall now examine the phenomenological consequences of the previous results for 
(double) layered superconductivity in parallel external magnetic field. Let us first consider 
the critical properties of the decoupling transition in the the absence of parallel magnetic 
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induction; i.e., take H\\ near H'^-^. Given the standard Ginzburg-Landau dependence, 
Al(T) = Ao(l — T/Tco)""*^/^, for the bulk {N — > oo) in-plane London penetration length, 
then (32) yields 

J||(T) = 7-i/cBTo(l-T/Teo) (56) 
for the intra-layer Josephson coupling energy, where 

is the basic energy scale of the problem. Since the zero-field decoupling transition occurs at 
/cfiT* = 47rJ||(T*), we then obtain T* = [T^^ -|- 7(47rTo)~-'^]~-'^ for the decoupling transition 
temperature. This means that the size of the critical regime is 

Sn = Tco -n= 77;o/47rTo (58) 

for Tco ^ To, which is typical. Likewise, the critical temperature at which each 
individual layer undergoes a superfluid KT transition is set by ksTc = f J||(Tc), or 
Tc = [T^^ + ^(^To)~^]~^. Hence, the distance of this intra-layer resistive transition to 
the Ginzburg-Landau transition temperature Tco is 5Tc = Tco — Tc^ 10 as indicated 
by Fig. 2. Again, we highlight the extraordinary regime in temperature < T < 
where the layers are normal yet Josephson coupled; i.e., py 7^ and H^}i{T) 7^ 0. This 
effect has been recently observed in resistance measurements on the highly anisotropic 
Bismuth-based series of high-temperature superconductors. ^^'^^ 

Results similar to those outlined above have been obtained recently by the author 
using an alternative anisotropic "frozen" superconductor model for the Meissner phase,® 
but with the important exception that the zero-temperature Josephson penetration length 
(8) appearing above in Eq. (57) is replaced therein by the zero-temperature London 
penetration length Aq. This discrepancy can be understood as follows: in the present 
frustrated XY model description (1) of the mixed phase, we take first the limit Aq — * 00, 
and then the limit ^ ^ 00, which results in the characteristic Josephson penetration 
length jd. In the anisotropic "frozen" superconductor model for the Meissner phase,® on 
the other hand, the order of the limits is reversed, hence the characteristic length scale Aq. 
Both models, however, obtain the same expression for the parallel lower critical field (33) 
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near criticality up to a numerical constant. In particular, we predict that H^-^iT) vanishes 
exponentially as temperature T approaches the decoupling transition from below, which 
implies the existence of an inflection point below in this temperature profile (see Fig. 
2)- 

In the presence of external parallel magnetic field, we expect that (double) lay- 
ered extreme type-II superconductors follow a cross-over phenomenon as a function of 
this field in the intermediate regime of the mixed phase. In particular, Eqs. (39) and 
(41) indicate (with Aq replaced by Aj) that the parallel magnetization has the form 
— 47rM|| = iyj,'j^(T)/[S||/s| (T)], where the latter functional dependence with parallel mag- 
netic induction is plotted in Fig. 1. Notice that in the mixed phase at low magnetic 
induction, B\\ <^ s|, we have that —A.ttM\\ ~ Hci{T). Hence, the parallel magnetization 
inherits the infiection of the parallel lower-critical field as a function of temperature T ~ 
at fixed S|| . Fig. 2 shows the phase diagram expected of a (double) layered superconductor 
in parallel external magnetic field near the critical regime discussed above. Formula (50) 
for the crossover field has been employed here, where the Josephson penetration length 
Aj(T) is interpolated between its behavior at criticality (31) and its low temperature value 
of 7(i. Notice that is expected to be practically constant at low temperatures since the 
mass anisotropy parameter 7 has no temperature dependence in this regime. We therefore 
expect the crossover field to exhibit an inflection point in its temperature profile, much like 
the parallel lower-critical fleld does. By Eq. (51), however, the ratio of to ijj^^ should 
be larger at criticality with respect to zero-temperature by a factor of A|(T*)/A|(0), which 
in the Ginzburg-Landau theory approximation is given by 47rTo/7T* ~ Tco/6T^. Last, in 
spite of the above crossover phenomenon, the parallel vortex lattice (and flux quantization) 
will persist up to the parallel upper-critical field. In the Landau- Ginzburg approximation, 
this field is set by the in-plane coherence length ^ — ^o(l — T /Tco)~^^'^ and by the mass 
anisotropy to be = 7^o/27r^^, hence the inequality bI <^ ffjg- The critical behavior 
of h12{T) near the decoupling transition at T*, however, remains unknown. 

We shall now examine the various physical scales that arise from the present the- 
ory in the context of high-temperature superconductivity. The oxide superconductor 
Bi2Sr2CaCu208 may be classified as a layered superconductor with an extreme mass 
anisotropy,-^ 7 ~ 100. Assuming typical parameters Tco ~ 100 K, d = 15 A and Aq ~ 10^ A, 
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we obtain from Eq. (58) that 5T^ ~ 0.5 K. It is interesting to remark that the zero- 
temperature Josephson penetration length and the zero-temperature London penetra- 
tion length, Ao, are of the same order of magnitude in this material. The fact that the 
present estimate for the size of the critical region is smaller by an order magnitude with 
respect to the estimate based on the above mentioned anisotropic "frozen" superconduc- 
tor model^ is simply then a result of the numerical factor in Eq. (58). Both theories, 
on the other hand, predict an inflection point in the temperature profile iyj,'j^(T) (see Fig. 
2). Wan et al. also observe an inflection point in the field of first penetration Hp vs. T 
for Bi2Sr2CaCu208, but in the perpendicular orientation.^^ This may be a vestige of the 
same prediction made here for the parallel lower-critical field if geometrical demagnetiza- 
tion effects are presumed to be strong. In particular, consider the regime in temperature 
Tc < T < T*, where the planes are resistive while remaining Josephson coupled, and hence 
where only parallel Josephson vortices exist. ^ Then the field of first penetration in the per- 
pendicular orientation is limited only by those portions of the field lines that run parallel to 
the top and bottom layers of the sample. Clearly, direct measurements of the parallel lower 
critical field in the critical regime of these materials would be highly desirable. Last, we 
mention that the parallel crossover field (50) at zero temperature, -Bi'(O) = {2/tt'^){^q/(P^), 
should be approximately 2 T for Bi2Sr2CaCu208, while it should be orders of magnitude 
smaller at temperatures just below the decoupling transition temperature. We therefore 
suggest that the parallel reversible magnetization be measured in a clean thin film of 
this material near criticality,^^ where the crossover field is expected to be quite modest 
in magnitude. Note that the present theory is valid only for perpendicular components 
of magnetic induction with magnitude greater than the perpendicular crossover field^^ 
~ ^o/d^l^i which in the case of Bi2Sr2CaCu208 is approximately 1 kG. 

V. Discussion 

In summary, we find no evidence for field dependence in the decoupling transition tem- 
perature, T*, of clean double-layered extreme type-II superconductors in the intermediate 
regime of the mixed phase. Since the double layer represents the weakest link, we believe 
that this result remains true in the general case of any finite number of layers as long as 
the interlayer coupling is weak enough so that the system remains effectively 2D; e.g., for 
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perpendicular fields above the Glazman-Koshelev 2D-3D cross-over scale, B± ~ ^q/^'^cP, 
which guarantees the absence of vortex loops that traverse many layers. -"^^ In general, the 
extreme type-II limit Nd <^ Xl must be taken first, however, so that magnetic screening 
efi'ects may be neglected. Also, the effect of pinning centers is not expected to be relevant 
in the critical regime, since the Josephson penetration length Aj(T) diverges exponentially 
as temperature T approaches from below. We do however find a cross-over field above 
which the parallel reversible magnetization decays with parallel magnetic field i^y as H^^ 
in the double-layer case (see Fig. 1 , Fig. 2, and ref. 27). The array of Josephson vortices 
is nevertheless expected to persist up to the parallel upper-critical field. It is important 
to mention that the present double-layer study cannot account for effects due to the first 
order commensuration transitions predicted to occur in the parallel vortex lattice with 
many layers. 

In order to account for the entropy due to wandering of the Josephson vortices in 
the calculations reported above, we have considered the length dimension of the vortex 
as imaginary time, and proceeded to compute the corresponding "quantum mechanical" 
correction to the "mass". Employing a renormalization-group improved semi-classical 
approximation to this end, we have found that the parallel lower-critical field vanishes 
exponentially as it approaches the decoupling transition temperature from below. Very 
similar results have been obtained recently by the author using an alternative anisotropic 
"frozen" superconductor model that operates from the Meissner phase. ^ Note that although 
the dimensionless Planck constant (17) has a value of hp — Sir at the decoupling transition, 
which is far from being small, it is suspected that the above cited renormalization-group 
improved semi-classical results are in fact exact for the case a single Josephson vortex. ^^'^^ 
Less is known, however, with respect to the validity of the semi-classical approximation in 
the critical regime for the case of the array. For example, we have computed the entropic 
pressure (49) of the array to first order in powers of the effective dimensionless Planck 
constant (17) and found it to be negligibly small at fields below a relatively large bare 
scale (43). However, we were unable to find a renormalization group for this contribution. 
Also, it has been argued in the literature that the pressure exerted between interfaces in 
two dimensions generally varies quadratically with temperature,^^ which translates into a 
second order correction in the present semi-classical approximation. These effects generally 
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stiffen the array of Josephson vortices and add a diamagnetic contribution to the high- 
field tail shown by the parallel magnetization (see Fig. 1 and ref. 27). However, any such 
additional diamagnetic correction does not affect the conclusion drawn here that no layer 
decoupling transition occurs as a function of external magnetic field in extreme type-II 
layered superconductors since it only results in a stiffer vortex lattice. 

The author is grateful to D. Dominguez, L. Bulaevskii, B. Ivlev, M. Maley, and H. 
Safar for very useful discussions. This work was performed under the auspices of the U.S. 
Department of Energy, and was supported in part by National Science Foundation grant 
DMR-9322427. 
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Figure Captions 



Fig. 1. Shown is the reversible parallel magnetization of an extreme type-II double-layer su- 
perconductor in the intermediate regime of the mixed phase; i.e., the parallel fields 
satisfy H^-^ <^ H\\ = B\\ -C h12- The cross-over field bI, in particular, satisfies the 
latter inequalities [see Eq. (51)]. Also, the tail in the magnetization that appears 
at fields beyond this scale varies asymptotically as Although these results are 

generally valid only for temperatures in the Ginzburg-Landau regime, they should 
provide a good lower bound for — 47rM|| at temperatures in the critical regime near 
(see ref. 27). 

Fig. 2. Shown is the phase diagram of a double-layer extreme type-II superconductor in par- 
allel external magnetic field near the critical regime. The KT transition temperature 
Tc marks the point above which each individual layer becomes resistive in the ab- 
sence of external magnetic field, while marks the layer-decoupling transition. Also, 
the dashed line represents the temperature profile for the parallel upper-critical field 
within the Ginzburg-Landau approximation. 
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